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Abstract: A continuum evolutionary model for microniagnetics is presented that, be- 
side the standard magnetic balance laws, includes thermo-magnetic coupling. To allow 

■^T ' conceptually efficient computer implementation, inspired by relaxation method of static 

• ' minimization problems, our model is mesoscopic in the sense that possible fine spatial 

,.^ I oscillations of the magnetization are modeled by means of Young measures. Existence 

of weak solutions is proved by backward Euler time discretization. 
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\Q • 1. Introduction, static problem and its relaxation 

. I Microniagnetics is a continuum theory introduced by Brown [S| to describe the equilibrium 

^l] • states of saturated ferromagnets. The equilibria are determined as minimizers of a functional 

with exchange, anisotropy, Zeeman (interaction), and magnetostatic energy contributions. This 
theory also predicts the formation of domain structures. The reader is referred to |17j for a recent 
survey on the topic. In the isothermal situation, the configuration of a rigid ferromagnetic body 
occupying a bounded domain fi C M'' is typically described by a magnetization m : Q ^ W^ 
S^ I which vanishes if the temperature 6 is above the so-called Curie temperature 6c and no external 

H ■ magnetic field is applied. 

- - -' On rriicroscopic level, the uiagnetic Gibbs energy consistsoiiivepa.rts,najne\y a.n anisotropy 

energy J^ip{m,9)dx, an exchange energy ^ J^e\\7m{x)\'^dx having a quantum-theoretical ori- 
gin, the non-local magnetostatic energy i /jgi/io|Vu,„(a;)pd2: = ^ J^^m-Vu„idx, an interaction 
energy — J^ h{x)-m{x) dx involving the outer magnetic field h and, finally, a calorimetric term 
/q ijjo dx. In the anisotropic energy, we denoted ij} its density that will depend on the material 
properties and should exhibit crystallographic symmetry. Furthermore, "tp is supposed to be a 
nonnegative function, even in its first variable. In the magnetostatic energy, m,„ is the magne- 
tostatic potential related to m by div(/xoVMm— Xo"t-) = arising from simplified magnetostatic 
Maxwell equations. Here xn '■ K'' — ?► {0, 1} denotes the characteristic function of il and /io is 
the permeability of vacuum. Moreover, it shall be noted that, due to the quantum mechanical 
origin, e in the exchange energy is very small from the macroscopic point of view. 

A widely accepted model for steady-state isothermal configurations is due to Landau and 
Lifshitz [20l[2T] (see e.g. Brown [5] or Hubert and Schafer [13|), relying on a minimum-of-Gibbs'- 
energy principle with as a parameter, i.e. 



minimize 



G^{m) :— I I '0(?Ti, 0) + — 771,-Vwm + 7:|Vto| ~h-mdx\dx 
Jo ^ 2 2 / 

subject to div(/ioVwm — XWrn) = in M^ , 



(1.1) 



where the anisotropy energy tp is considered in the form 

4>im, e) :- 0(m) + ao{e - e,)\m\'' - M^), (1-2) 

where gq determines the intensity of the thermo-magnetic-coupling. To see a paramagnetic state 
above Curie temperature 0c, one should consider oq > 0. The isothermal part of the anisotropy 
energy density (j> : R'' — > [0, cxd) typically consists of two components (j){m) = 0poios("^) + ^o|'7i|'*, 
where (/'poies('7i) is chosen in such a way to attain its minimum value (typically zero) precisely 
on lines {tSa, t € M}, where each Sq, G M.'^, \sa\ = 1 determines an axis of easy magnetization. 
Typical examples are a = 1 for uni-axial magnets and l<a<3orl<a<4 for cubic 
magnets. On the other hand, 6o|m|^ is used to assure that, for 6 < 9^, V'('7^) is minimized 
at tsa for |t|^ = {9c — 9)a^/{2ho). Such energy has already been used in [31]. For e > 0, the 
exchange energy ejVwp guarantees that the problem (jl.ip has a solution m^. Zero-temperature 
limits of this model consider, in addition, that the minimizcrs to (jl.l[) are constrained to be 
valued on the sphere with the radius y^ao9c/{2bo) and were investigated, e.g., by Choksi and 
Kohn [7| , DeSimone [5] , James and Kinderlehrer [T3| , James and Mliller [TS] , Pedregal [551 112] : 
Pedregal and Yan [30] and many others. 

For £ small, minimizers rUe of (jl.ip typically exhibit fast spatial oscillations, a so-called 
fine structure. Indeed, the anisotropy energy, which forces magnetization vectors to be aligned 
with the easy axis (axes), competes with the magnetostatic energy preferring divergence-free 
magnetization fields. If the exchange energy term is neglected, and this is a justified simpli- 
fication of the functional for large ferromagnets [H], nonexistence of a minimum for uniaxial 
ferromagnets can be expected and was shown in [T3] for the zero external field h and zero 
temperature. Hence, various concepts of relaxation (in the sense of variational calculus) were 
introduced in order to cope with this phenomenon. The idea is to capture the limiting behavior 
of minimizing sequences of G^{m) as e — > 0. This leads to a relaxed problem (|1.3|) involving 
so-called Young measures i^'s which describe the relevant "mesoscopical" character of the fine 
structure of m. We call this "limit" a microstructure. 

It can be proved [^I28j that this limit configuration (v, u„i) solves the following miniiniza- 
tion problem involving "mesoscopical" Gibbs' energy G: 

minimize G{v, m) := / yip ti/ + —m-'Vuni — h-iri) Ax 
subject to div(/^oVum - xo^) = on M'', i , , 

TO = id»i^ on il, 

i/£^P(^;R'*), meLP{Vl]W^), u,neH\R'^) , 

where the "momentum" operator " • " is defined by [f *i']{x) := J^^ f{s)i'x{<is), id : R"* — s> R"* 
denotes the identity and v G ^P(fi;R'*). Here, the set of Young measures ^P(r2;K'^) C 
L^{fl]M{R''')) = L^i^; Co(R'^))* is the set of aU weakly measurable essentially bounded map- 
pings X ^ Vx : n ^ 7W(R'') = Co(R'')* such that J^J^^ \s\Piyx{ds)dx < +oo; here Co(R'^) 
denotes the set of continuous functions with compact support and thus A^(R'') is the set of 
Radon measures on R'', and the adjective "weakly measurable" means that vi^ is Lebesgue 

measurable for any v e Co(R''). A natural embedding i : LP(fi;R'*) -^ ^P(fi;R'') of a mag- 
netization 771 is a Young measure v = i{m) defined by i>x = Sm{x) with Ss denoting the Dirac 
measure at s 6 R''. We say that a sequence {:^'^}fcgN C ^P{il;R.'^) converges weakly* to v if 

V/eLi(0;Co(R'')) : lim ff.iy''dx= ff.i^dx (1.4) 

The set ^P(r2;R'*) is convex, metrizable, compact, and contains weak*-densely the set of mag- 
netizations TO G L^(i7;R'') if embedded via i. One can thus call '3^p{Q,]W^) a convex compact- 
ification of the set of admissible magnetizations, cf. also j331 Chapter 3]. The so-called relaxed 
problem p.3p can then be understood simply as a continuous extension of the original problem 



(|1.1|) considered for e = 0. Let us note that the problem (|1.3p has a convex structure. Moreover, 
it captures the multiscale character of the problem. The Young measure solving (|1.3|) encodes 
limiting oscillating behavior of minimizing sequences of (|1.1|) while its first moment, the magne- 
tization, resolves the macroscopic magnetization m. Due to these properties we call the model 
in (|1.3p mesoscopic. An equivalent way how to relax (|l.ip is to replace tJj by its convex enve- 
lope. The drawback for numerical calculations is that one needs to know the convex envelope 
explicitly. There were many attempts to design numerical schemes for both (jl.ip , as well as for 
(|1.3p in the zero-temperature situation; cf. e.g. [6l[l6l[22]- Departing from p.3p and following 
the ideas of [26], a model of the isothermal rate-independent evolution exhibiting hysteretic 
response was proposed and analyzed in [35, 36J. 

Our motivation is to merge concepts of relaxation that can successfully fight with multi- 
scale character of the problem with recent ideas to build thermodynamically consistent meso- 
scopic models in anisothermal situations. A closely related thermodynamically consistent model 
on the microscopic level was introduced in j31j to ferro/para magnetic transition. Another re- 
lated microscopic model with a prescribed temperature field was investigated in [1] . The goal is 
to develop a model that would be supported by rigorous analysis and would allow for compu- 
tationally efficient numerical implementation like in |16| I18| I19| where such a model was used 
in the isothermal variant. 

This paper is structured as follows: In Section 2 we introduce the concept of general 
standard materials and handle thermodynamics, in Section 3 we give a weak formulation of the 
model equations, in Section 4 we give a time-discretization of the model equations and prove 
a priori estimates, in Section 5 we then prove convergence of the discrete solutions. Finally in 
Section 6, we present a generalization of the model, where we allow for a weaker dissipation so 
that at least some internal parameters can be considered evolving in a purely rate-independent 
manner. 



2. Evolution problem, dissipation, mesoscopic Gibbs free energy 

If the external magnetic field h varies during a time interval [0,r], T > 0, the energy of the 
system as well as the magnetization evolve, too. Change of magnetization may cause energy 
dissipation. As the magnetization is the first moment of the Young measure we relate the 
dissipation on the mesoscopic level to temporal changes of some moments of v and consider 
these moments as separate variables. This approach was already used in micromagnetics in 
[35[l36j and proved to be useful also in modeling of dissipation in shape memory materials, see 
e.g. [24]. In view of ()1.2p . we restrict ourselves to the first two moments defining A = (Ai, A2) 
giving rise to the constraint 

\ = L»v, where L(m) := (m, |mp) (2.1) 

and consider the specific dissipation potential 

C(A):=<5J(A) + J|Ar', q>2, (2.2) 

where Sg : M.'^'^^ — ^ R"*" is the Legendre-Fenchel conjugate function to the indicator function 
5s ■ M.'^'^^ — !■ {0, +00} to a convex bounded neighborhood of 0. The set S determines activa- 
tion threshold for the evolution of A. The function Sg > is convex and degree- 1 positively 
homogeneous with Sg{Q) = 0. In fact, the first term describes purely hysteretic losses, which 
are rate-independent (the contribution of which we consider dominant) and the second term 
models rate-dependent dissipation. 



In view of (|1.2P " (|1.3p . the specific mesoscopic Gibbs free energy, expressed in terms of u, 
X and 0, reads as 

g{t, iy,X,9) ■.^(j)*i^ + {0-0c)a-\ - tpoiO) + -m-/idcm - h{t)-m (2.3a) 

with m = id»h' and /idem = Vum, (2.3b) 

where we denoted a := (0, . . . , 0, uq) with aq from ()1.2|) and, of course, Um again from p. ip .which 
makes g non-local. 

In what follows, we relax the constraint (|2.ip by augmenting the total Gibbs free energy 
(i.e. g integrated over fi) by the term -f ||A — i • i^||^_i,j^.jgrf+i, with (presumably large) >c G M+ 

and with H^^{rt; ]R'^+^) = ^o (^; 1^"^+^)*. Thus, A's no longer exactly represent the "macroscop- 
ical" momenta of the magnetization but rather are in a position of a phase field. Let us choose 
a specific norm on H-^{n;R'^+^) as ||/||ff-i(n;R'i+i) := \\x\\L2{n;R''+^) with x G -^gradl^;'^'') 
defined uniquely by divx = X-L*v where L^^^^C^;'*'') — {^«; " ^ i/o^(fJ;K'^+i)}. In other 
words, ||/||H-i{f2;i{d+i) = II VA~"'^/||^2(f2.jjd+i-) with A^-' meaning the inverse of the Dirichlet 
boundary-value problem for the Laplacean A : _ffQ(17;IR^+^) — > H^^{fl;R'^^^) . We define the 
mesoscopic Gibbs free energy '^ by 

'^{t,,^,X,0) :=y (5(t,z^,A,0) + ||VA-i(A- L.j.)p) dx (2.4) 

Notice that the i?~^norm defined above is equivalent to the standard i?^^-norm defined 

as ||/||_i ;= sup,^,g^y^l,2(J^.JJ^+l^_ ^_^^ {f,v) /\\v\\^yi.2^^.^^+,y Indeed, if /i £ Wo'^{n;W^+^) solves 
-Ah = f and \\v\\yyi.2^^.^i+i^ ~ ||Vu||L2(o.Rd+i) then 

11/11-1= sup ll«ll"\.2,„„^+i. / V/t- Vi>d.T >/3||V/i||i2(n;Rd+i) , 

where /3 > is the ellipticity constant of B[h, v] := J^^ \7h ■ Vv dx. On the other hand, 

11/11-1= sup ||Hlwl^ro.».+n /v/^-V«dx<||V/i|U 



dSWo' (n;R<'+l), v^O 



wi-^-{n;m''+^) / ^"' ^''"^' - ll^"lli^(n;R''+i) 



by the Cauchy-Schwartz inequality. 

We shall give a certain justification of the penalized model at the end of the section. Yet, 
we should emphasize that we will consider >c fixed thorough this article. 

The value of the internal parameter may influence the magnetization of the system (and 
vice versa) and, on the other hand, dissipated energy may influence the temperature of the 
system, which, in turn, may affect the internal parameter. In order to capture these effects 
we employ the concept of general standard materials [12| known from continuum mechanics 
and couple our micromagnetic model with the entropy balance with the rate of dissipation on 
the right-hand side; cf. (|2.6p . Then i/ is considered evolving purely quasistatically according to 
the minimization principle of the Gibbs energy Q{t, •, A, 9) while the "dissipative" variable A is 
considered as governed by the flow rule in the form: 

aC(A) = 9A5(i,^,A,0) (2.5) 

with 9C denoting the subdifferential of the convex functional (■(■) and similarly d\g is the 
subdifferential of the convex functional g{t,v, -,9). In our specific choice, (|2.5p takes the form 

d5g{X) + e\X\'^^^ X+ {9—9c)a 3 m:A^^{X — Ltiy). Furthermore, wc define the specific entropy s by 

the standard Gibbs relation for entropy, i.e. s ~ -~g'g{t, v, A, 0), and write the entropy equation 

9s + div j = ^(A) = heat production rate, (2.6) 



where j is the heat flux. In view of (|2.2p , 

e(A)-aC(A).A-5*(A) + 6|Ar. 
Wc assume j is the heat flux assumed governed by the Fourier law 

j = -KV6' 



(2.7) 



(2.8) 



with a heat- conductivity tensor IK = IK(A, 9). Now, since s = ~g'g(t, v, A, 9) = —g'g{X, 9), it holds 

9s = —9gg{X, 9)9 — 9gg^\. Using also g'g^ = a, we may reformulate the entropy equation (|2.6p 
as the heat equation 



c^{9)9-diw{K{\,9)V9)^5*s{\) + e\X\'' + d-9\ with c^{9) ^ -9g'^{9), 

where Cv is the specific heat capacity. 

Altogether, we can formulate our problem as 



(2.9) 



minimize 



(4>^v+{9(t))^9,)a-\{t)-i,o{9(t)) + -m-/idcm 



-/i(t)-m +f|VA-i(A(t)-i. 



d.T 



on O, 



► for ie[0,T], 



subject to m ~ \(i*v, /idem = Vu„i 

div(/^oVum - Xfi'™) = on M'', 

z^G^P(f7;M''), m(zLP{VL-W^), u^eH^iR'^), 
dSsiX) + e| A|«-2 A + {9-9,)a 3 hI\-^ (div X-L.v) in Q := [0, T] xfJ, 

c^{9)9 - div(K(A, 9)V9) = (55(A) + e|A|« + a-A6i in Q, 

(K(A, 6*) V6I) -n + b9 = b9^^t on E := [0, T] xT, 



(2.10a) 



(2.10b) 

(2.10c) 
(2.10d) 



where we accompanied the heat equation (|2.10cp by the Robin-type boundary conditions with n 
denoting the outward unit normal to the boundary F, and with b G L°°{T) a phenomenological 
heat-transfer coefficient and 9cxt an external temperature, both assumed non-negative. 

Next we shall transform (j2.10c[) by a so-called enthalpy transformation, which simplifies 
the analysis below. For this, let us introduce a new variable w, called enthalpy, by 



w ~Cy{9) = / Cv(r)dr, 
Jo 



(2.11) 



It is natural to assume Cv positive, hence Cv is, for w > increasing and thus invertible. 
Therefore, denote 



eiw) 



Cy ^ (w) if u' > 
if w < 



and note that, in the physically relevant case when 9 > 0, 9 ~ Q(w). Thus writing the heat flux 
in terms of w gives 

K(X,9W9^K(X,e{w))VQ(w)^JC(X,wWw where IC(X,w) := ,'^,2' . (2.12) 

Cv(B(w)) 

Moreover, the terms {Q{w{t))—9c)d-X{t) and tpo{9{t)) obviously do not play any role in the 
minimization (|2.10ap and can be omitted. Thus we may rewrite (|2.10ap in terms of w as follows: 



h(t)-m+^\VA-\X{t)-L.u)\^)dx 



mininiize 

subject to m ~iA*v, /idem = "^Um on i7, 

div(^oVM,„ - Xf2™) =0 on K'', 

9<5J(A) + e|A|«"2A + (e(«;)-e'c)a 3 kI\-~^{\^L,v) in Q, 

w - div(/C(A, u;)Vu;) == (55(A) + e|A|« + a • GHA in Q, 

(/C(A,w)Vw)-n + 6eH ^^^eioxt on E. 

Eventually, we complete this transformed system by the initial conditions 

z^(0,-)=i'o, A(0,-) = Ao, ^(O, •) = Wo := Cv(6'o) on Vt 



>forie[0,r],(2.13a) 



(2.13b) 

(2.13c) 
(2.13d) 



(2.14) 



where {vq,Xq) is the initial microstructure assumed to solve (|2.13ap and the phase field, and 
00 is the initial temperature. Note also that, by prescribing j/q, we also prescribe the initial 
magnetization toq = id • j/q and magnetic potential Umo ■ 

2.1. Justification of the penalization concept 

Recall that in the model (|2.10p we gave up the constraint \ = L*v and only included a 
penalization term ^\\\ — i»i^||_y-i(f2:K<'+i) in the Gibbs free energy. To justify this approach, 
we show that in some particular situations, namely in the static case and also in the iso-thermal 
rate-independent (with a small modification) case, solutions of the penalized model converge 
to solutions of the original model that satisfy A = L*v as x ^ oo. We shall also give some 
heuristic ideas, why a similar limit passage should be possible even in p.lOp . however a rigorous 
proof is beyond the scope of this paper. 
The static case: 

Let us consider an analogical problem to (|1.3p that includes also a penalization term and 
where we also use the form of the Gibbs free energy as in (|2.3ap with 9 given, i.e. 

1 X 

0.1/ + {9-0c)d-X - V^o(^) + -m-Vw,„~ft,(i)-m+— |VA"^(A - L.j/)p) d: 



<2.15) 



minimize 

subject to div(/xoVum — Xn^ ~ on M , 
TO = id»z^ on J7, 

iye'3^P{n;R'^), A G i/-i(0;M'*+i), meLP(r2;]R''), u,„ei/i(K'') 

Let us denote {X^,iy^) G H-^{n;W'+^) x '3^P{n;W^) the solutions to ([235]) . Let us then 
show that they (in terms of a subsequence) converge weakly* in H~^{il;M.'^^^) x ^P(r2;IR'^) to 
the solutions of 

1 



minimize / l(f>*i'+{0—9c)d-X — ipo{0)-\ — m-'Vum — h{t)-m] dx 
subject to div(/ioVum — xn^ =0 on M , 
m~id*i/, X — Lmv on i7, 

u€'3^P{n;R'^), Xe H-^{n;M.'^+^), TOeLP(f7;R''), w„,Gffi(M'')J 



>(2.16) 



Namely, it easy to see, from coercivity of (j) and by simply testing p.lSp by any (A, O) such 
that X = L*v, that /^^ | • |^ • v^^dx is bounded uniformly with respect to k; here p corresponds 
to the growth of 0, cf. p. 5b ,). Hence also ||i»i/><||L2(n:R''+i) and in turn also ||L»z/>f||/f-i(o.R<i+i) 
are uniformly bounded with respect to x. 



Also, by the same test as above, we get that >c\\\^ — L»v^\\^_i,^,^i+i-. and thus also 
||A>c||^_l(•Q.J^d+l^ are bounded independently of x. 

Hence, exploiting standard selection principles, we find a pair {v^\) G ^^(r2;M'') x 
H~^{Vl\W^^^) such that (in terms of a not-relabeled subsequence) i^^^ ^ i^ in '3^P{Vt]W^) and 
Aj< ^ A in H~^{VL\W'-^^). Also, as k\\\^ — i»J^><||^-i(Q.Rd+i-| is bounded independently of x, 
necessarily X = L»v holds for the weak limits. 

Then thanks to the weak-lower semi-continuity we have that 

{4)»v+{9—9c)a-\—tpo{9) + —m-yum—h{t)-rn)j dx 

< liminf / [(f)*v^-\-{9—6c)a-\^~'tl}Q{9)'{'—m^-'SIum —h{t)-m^) dx 

< liminf / (<p»i>+{9-9c)d-X~ipQ{9) + -rhVurh~h(t)-7fi] dx+—\\\-Lti>\\jr-i(a.Bd+i) 
= / ((j)^i)+{9—9c)a■X—^|Jo{9) + -rhVu,n—h{t)■7fljdx 

for any (A, D) such that \^ L*v, which shows that (A, v) is a solution to (j2.16p . 

The rate- independent isothermal case: 

When considering the rate- independent case we formally set e = in (|2.2p . i.e. assume 
that the dissipation potential is equal to 5*g. Also, since now the dissipation potential yields 
less regularity on A we have to alter the specific mesoscopic Gibbs free energy and to add a 
regularization term 7|VAp; let us therefore denote 

gRi(t,i.,A)=.g(t,i.,A,0)+7|VA|2, 

again for some 9 fixed. Now, as the temperature is not a variable in this context, the system of 
governing equations (|2.10p reduces to (with 5Ri(i, v. A) replacing g{t, i/. A, 9)) 



minimize (gB_i{t,i',X,)^ IVA ^{\{t)—Lti')\ ) dx 

subject to m~[d*i>, on 51, 



>ioi-te[0,T], (2.17a) 



div(^oVu„ - Xnm) =0 on M , 

ve?VP{n;R'^), meLP{n-R'^), u,„eHi(K'^) 
dS*g{X) + {9-9c)d + 2-/div\7 X + xA-^ {div {\ - L.i^)) 3 in Q := [0,T]xSl, (2.17b) 



and we recover, apart from the penalization, a similar model as in |18 1 I19 ) . Since we are consid- 
ering the rate- independent case, a suitable weak formulation of (|2.17p . which in fact - due to 
the convexity of the problem - is equivalent to the standard weak formulation, is the so-called 
energetic formulation, cf. e.g. [26]. Then we shall cah (A^,i/^) G L°°{[0,T];W^^'^{n;R'^+'^)) n 
BV([0,T];ii(r2;R''+i)) x (^P(fi;ffi''))[o^'^l an energetic solution of (^37]) if they satisfy (we 



included initial conditions here already) 

< f (gH,i(i,:), A) + ||VA-i(A-L.i>)|' + <5J(A - A^))dx 

for all (A, i>) G W^^^{n- R'^+^) x ^^(17; W^) and all t G [0, T] , (2.18a) 

(.gRi(r,;.^(T),A^(T)) + ||VA-i(A^(r)-L.i.^(r))|')d.T + Var55(A;0,r) 



< / .gRi(0,i^o,Ao)dx+ / / [gRi];(s,i^^(.s))dxds, (2.18b) 

Jn Jo Jsi 

X^{0)=XoeW^^\n;M.'^+^), iy^{0) = uoe3^P{n;R'^), Xo = L,i^q, (2.18c) 



with Var/(x;0,T) the space integral of the variation of / between and T. Let us now 
show that energetic solutions (A^,z^^) of ([STT]) converge in L°°{[0,T];W^'^{n;R'^+^)) fl 
BV{[0,T]; L\n;R'^+^)) x (^P(fi;R'*))[o^'^l (at least in terms of a subsequence) for >f ^ cx) 
to (A, v) satisfying 



gm{t,iyit),X{t))dx< / (gRi(i,i),A) + 5S(A-A))dx 
for ah (A, i>) G W^'^{n; R'^+^) x ^^(1]; R'^) such that A == L . i> and alH G [0, T] , (2.19a) 



(7ai(T,i/(r),A(r))da; + Var5j(A;0,T)< / g^j{0,,yo,Xo)dx + / [g^j]'t{s,i^{s))dxds, 
n Jn JQ 

(2.19b) 

A(0) = Ao G VF^'2(f7;R''+i), i^(0)=z^o e^P(f};lR''), Ao^L.i^o, (2.19c) 

i.e. the energetic formulation of 



minimize / gKi{t, v^ A) dx 
Jn 



in 
subject to ?Ti = id • i^, X = L »v, on Q, 



>foriG[0,r], (2.20a) 



div(^oVwm - xnm) =0 on M'', 

ve'3^P{n;R'^), m£LP{n;R'^), u^GiJ^M''), J 
dds{X) + {0-ec)a + 2-f^\7X + dIx=L.u3O in Q := [0,r]xfl, (2.20b) 

where Ix=l mu = 0, if A = L • i^, and +oo otherwise. 

The conjecture follows from the abstract paper |25| : here we give a very short sketch. 
By similar tests as in the static case, it can be seen that 

l|Ax:||L~([0,T];Wi.2(n;R<i+i))nBV([0,T]:Li(O;R'i+i))j Jn\'\^ ' ^'<'^^ ^'^"^-^ >^||A,j ^ ^ ' ^^\\H-^{n;R'i+^) ^'"'^ 

bounded independently of k. 

Hence, by a slight modification of HcUy's theorem |231 [25] . there exists A G 
i°°([0,T];M^i'2(rj;M''+i)) nBV([0,T];ii(f7;M''+i) and a not-relabeled subsequence of x such 
that A^(i) -^ X{t) for all t G [0, T] weakly in W^^^{n; R'^+^) and hence strongly in L'^{n; R'^+^)). 

Let us now fix i G [0, T]. Then there exists a subsequence of x (dependent on t) denoted 
Kt and a 1/ G 3^P{fl;R'^)) such that z^^, -^ v weakly* in ^P(fi;R''). Similarly as in the static 



case, using (j2.18ap and weak lower semi-continuity, we get that 

gni(t,v{t),\{t))dx < liminf / gmit,v^tit), X^,it))dx + ^\\X^^ - ^ • '^>^JlH-i(n-R<i+i) 

< liminf y (gRi(i,i>, A) + S^'^ - A^(t)(i)))dx = J (.9Ri(i,i>, A) + S*sCX - A(i)))dx, 

for all (A, v) e W^'^^fl; R'^+i) x ^P(rj; R'^) such that \ = L.v, i.e. we showed that (A, v) fulfills 
()2.19ap . In the last line we exploited that 5%, is one-homogeneous and the strong convergence 
\^{t)^\{t) {nL^{n-W^+^)). 

To see that (A, v) also fulfills (j2.19bp we pass to the limit in (j2.18b[) exploiting only weak 
lower semi-continuity. 
The thermally coupled case as exposed in (|2.13p ." 

As already mentioned, we only give a short heuristic sketch why the penalty approach is 
also justified in the case presented here, in particular we concentrate only on the limit passage 
in the minimization principle p.l3ap since this seems to be the most involved one. Let us for 
simplicity assume that q ~ 2. 

Assume that (A^,i^^) e W^'^{[0,T];L^{n;R'^+^)) x (^p(J7;M''))[o^'^1 (together with some 
w^ S L^{[0,T];W'^'^{il)), which is however irrelevant here) are weak solutions of (|2.13p with 
an initial condition satisfying Aq = L* vq. Then (|2.13ap yields, just by the chain rule, that 

/ ((/).i'^(i) + ^™^W-Vu„,^(i) - h{t)-m^{t) + ||VA-i(A^(i)-L.zy^)|')da; 

= / ( (p' 1^0 + -;^'rnQ{t)-Vu,no{t) - h{t)-mo{t)]dx 

{\7A-\X^{s)-L.iy^{s)))-{VA-^X^{s)) - /i(s)-m(s))da;ds 



/o Jn 



for any t G [0,T]. Combining this with the fiow rule tested by A>c gives that supj^rp t] Jn\ ' 
\^ 'V^ is bounded independently of x and, moreover, that X^ is bounded independently of x 
in M^^'^([0,T]; L^(f2;M''+^)), too. Using the estimates for X,^ once again in the flow-rule, we 
get that /p^ |>!:/j-j VA"i(A^ - i . v^)VA^^vdx\dt is bounded for aU v S L'^{Q] W^+^) such that 
lkllL2(Q;K<i+i) < 1, in particular {\x j^WlS.^'^(Xtt{t) — L»i':^{t))'S/lS.~'^VQ,dx\}^yQ is bounded for 
all Us G i^(f2;M''+^) such that ||ws||L2(n.Rd+i') < 1 and a. a. t G [0,T]. This in turn means that 
{>f||A>f(t) — L * v^{t)\\H-i [ns.'i+i)} >c>Q is bounded for almost all t G [0,r]. 

Using that H/i^2([o, T];L^{n; W^+^)) C C([0, T]]L^{n; W^+^) one can select a subsequence 
of X (not relabeled) and find A G W^^'^([{),T];L'^{D.]W^+^)), such that X^{t) -^ X(t) strongly in 
H~\n;W^+^) for alHG [0,T]. 

Let us fix some t (1 [0,T] such that {K\\X^{t) ~ L ^h'^{t)\\ff-l|^Q.Jg^d+l^}:^ryo is bounded (note 
that this is possible a. a. t G O). Then similarly as in the isothermal rate-independent case one 
can find a subsequence of k (dependent on t) denoted xt and i/ G ^^(17; M'') such that ly^^ — > 
1/ weakly* in ^^(J7;M''). Once again by lower semi-continuity one can get the minimization 
principle 

/ ( i/)*;^— 7Ti-Vu,„— /i(t)-77i J dx < liminf / {cf>»i'3^t +—m^_^-Vu„i ~h{t)-m^^)dx 
Jo^ 2 / xt->oo 7j^ V 2 "' 

< liminf / ( 0. i/^,+-m^,-Vu„ -/i(t)-m^J da;+-^||A^,-L . i^^JI^-w^^.g^+i) 

^t^ao Jj-j V 2 2 y ' I 

< liminf / ((p*i'+-mWufn-h{t)-rh)dx+—\\X^Jt)-L»i)\\jj-i,^.^d+i) (2-21) 
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for any v e ^''(i7; M'^). It would seem logical to take v such that A(t) ~ L • D, then the 
penalization term on the right-hand side would become >c\\Xj^(t) — ^{t)\\'jj-i(Q.-jg^d+i)- However, 
although we know that ||A>c(i) — X(t)\\'ir-i (n-md+i) converges to as >f — > cxd this no longer needs 
to hold if the term is multiplied by >f. 

The limit passage in the minimization principle in the thermally coupled rate-dependent 
case therefore seems to be much more involved than in the cases presented before. This is mainly 
due to the fact that now the evolution of A is given by a completely separate equation. However, 
we can use a trick to circumpass this problem. Namely we realize that, due to the convexity^ 
any solution of (|2.13ap ly^, solves also the following problem (for the fixed t) 

minimize / (<fyv+ -m-Vu^ - h{t)-m + >i:{V A^^ {\^{t)-L m ly^)) ■ (VA"^(A^(i)-L . i^)) j dx 

subject to m = idmi/, /idem ~ Vum i^K solution to (|2.13ap on fi, 

div(/^oVMm - Xnm) =0 on M'', 

^e^P(f7;R'*), meLP{n;R'^), u^eH^iR"^). 

(2.22) 
Now we return to the second line ()2.2ip (and use that >ff /2 < >ct) but instead of exploiting 
that i/^^ solves ()2.13a|) we use that it solves p.22|) and get 

liminf / ((l)»iy^^ + -m^^-Vu,n-h{t)-m^^)dx+}<t\\\^^-L*iy^^\\]j-i(^.^d+i) 

= lim inf / ( . i^^, +-m^j ■ Vu„ ~h{t)-m^^ -h><t V A"^ (A^, -L « i^^J VA"^ (A^, -L »v^M dx 

< liminf / ( 0»£'+— mVuTO~^(i)''7i-t->ftVA~ (Aj^^— L»i/j^jVA^ (\^^—L»ii)\dx 

for all ii G ^P(ri;K'') such that A = L*v. Now we exploit our special choice of t for which 
we know that >ft||Aj^j(i) — i»!^3<t(i)||_ff-i(a;R<'+i) is bounded and ^>ttk'^) -^ A(t) strongly in 
_ff^^(i7;R'^+^), which makes the penalization term vanish. Hence we establish the converged 
minimization principle. 



3. Weak formulation, qualification and main results 

In this section we shall give a weak formulation of the proposed micromagnetics model. The 
used formulation is to a great extend inspired by the energetic formulation for rate- independent 
processes (see e.g. [IHl) and its generahzation given in e.g. [31] for problems that include both 
rate-independent and rate-dependent processes. 

In our case as well, we may regard the magnetic variable i> to be fast evolving; its evolution 
is therefore driven by rate-independent processes. On the other hand, the variables and A are 
evolving slowly as there evolution is driven by a rate-dependent process. Hence, we demand 
the fast evolving variable to satisfy a minimization principle and for the slow variables we just 
require standard weak formulation of (|2.13[ b-d) . 

Definition 3.1 (Weak solution). The triple {v, A, w)£{'3/'p{VI- M''))[0-^1 xiyi'«([0, T]; L«(r2; W^+^)) 
xLi([0,r];Wi'i(f})) such that m^ id* v <E L'^{Q]W^) andL.ue L'^{Q;R'^+^) is called a weak 
solution to (|2.13p if it satisfies: 

1. The minimization principle; For all v in ^P(ri;R'^) and all t G [0,T] 

^(t, V, A, Q{w)) < ^(t, i>. A, e{w)). (3.1) 
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2. The reduced MELxwell system for magnetostatics; For a. a. i G [0, T] and all v G H^{M. ) 



fiQ / Vum'Vwdx= / TO-Vwdx. (3-2) 

3. The flow rule; For any v e L9(g;M''+i) 

f neH-6'c)a-(w-A)+<5S(w) + -|t;|« + >fVA-i(A-L.i^)-VA-i(v-A)j dccdt 

> / (S*si\) + -IX]"^] dxdt. (3.3) 

4. r/ie enthalpy equation; For any ip G C^{Q), (p{T) = 
I /C(A, w)Vit;-Viy9 — ui(^ j dxdi + / bld{w)(pdSdt ~ / it;o<y2(0) da; 

('5S(A)+e|A|« + e(w)a-AVda;dt+ / bO^^tfdSdt. (3.4) 



s 
5. The remaining initial conditions in (|2.14p ; 1^(0, ■) = i^o and A(0, •) = Aq. 

Data qualifications: 

Let us now summarize the data qualification, needed to prove the existence of weak solutions: 

Isothermal part of the anisotropy energy: (j) G C(R'') and 

3c^c^>0, p>4: ct{l + \-n<H-)<4a + \-n, (3.5a) 

Rate-independent dissipation: 6g G C{M.'^^^) positively homogeneous, and 

3ci,z3, C2,i5 > : ci,z5(| • I) < 5*s{-) < C2M\ ' I), (3-5b) 

External magnetic field: 

heC\[0,T];L^{n;m.'^)), (3.5c) 

Specific heat capacity: Cv G C(R) and, with q from (|2.2p . 

3ci,e,C2,e > 0, cji > w > g', ci,e(l+0)""' < Cv(0) < C2,e(l+e)"^"', (3.5d) 

Heat conduction tensor: /C G C{W^+^ x M;R''^'') and 

3Ck,a.o>OVxgM'^:/C(-,-)<Ck, X^/C(-, Ox > «o|xP, (3.5e) 

External temperature: 

^cxtei'(S), ^oxt>0, and 6Gi°"(I]), 6 > 0, (3.5f) 

Initial conditions: 

1^0 e^''(^;K'') solving dSIIl) , Ao G L«(0;R''+i), iuq = Cv(6'o) e L^(fi) with 6^0 > 0. 

(3.5g) 

Note that (|2.1ip combined with (|3.5dp yields for non-negative w 



w=Me)= / Cv(r)dr>ci,e / (l+r)"-Mr > ci,(,((l+0)"- l) = ci,(,((l+e(u;))"- l). (3.6) 
Jo Jo 

The main analytical result we will prove in the following sections is: 
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Theorem 3.2. Let (j3.5p hold. Then at least one weak solution {ly, A, w) to the problem (|2.13|) in 
accord with Definition \3.1\ does exist. Moreover, some of these solutions satisfies also 



weL'-{[0,T];W''''{n)) n W'''{I;W''°°in)*) withl <r< 



d+2 
d+1' 



(3.7) 



Scenario of the proof. The above assertion will immediately follow from the Proposition 14.31 
which assures convergence of the approximate solutions constructed by semi-implicit time dis- 
cretization in Section 2] to the weak solutions of (|2.13|) , when proving a-priori estimates in 
Proposition 14.21 and when realizing that the approximations wq.t and Aq.t of wo and Aq re- 



quired in (|4.2ap below always exist. The information ^ g L^ {I;W^'°° (fl)*) can be obtained 
from the equation (|2.13b .d) itself. D 



4. Time discretization, a-priori estimates, and convergence 



To prove Theorem 13. 2[ wc proceed in a constructive manner that may serve also as a concep- 
tual numerical algorithm, at least after a spatial discretization being performed. Namely, we 
discretize the time with a time-step r and introduce a minimization problem in every time- 
step that we shall call the time-incremental minimization problem. This problem represents 
a discrete version of the minimization principle as well as of the flow rule. Also, we apply a 
semi-implicit method of time-discretization in such a way that it decouples the time incremen- 
tal minimization problem and the heat equation in any particular time-step k by using the 
"retarded" enthalpy, i.e. w^~^. 

We call the triple (^'^,A^,^y^) € ^P{n;R'^) x L^'?( »;R^+i) x H\n) the discrete weak 
solution of (|2.13p subject to boundary condition (|2.13dp at time-level k, k = 1...T/t, if it 
satisfies: 

1. The time-incemental minimization problem with given A^^^ and w'^"^: 

'A-Ai;-i\ e A-A^-^ 



Minimize '^{kT,iy, X,e{w';-^)) + t |Ap9 + 5 



dx 



(4.1a) 



subject to (t/,A) e '3^P{n;W^) x L^'i{n;W^+^). 



with^ from (plij) . 
2. The reduced Maxwell system for magnetostatics: For all v G H^(M.'^) 



yUmk -Vf dx 



m-Vvdx 



with m2; = id • ly^ 



(4.1b) 



3. The enthalpy equation: For all ip G H^{n) 






(p + /C(A^,<)V<-V^l da;+ / 6^e(<)^dS' 
x'i-x':-^\ Ai:-A^-i 



fe^^t.^y'ds 



's 



e(w^)a 



4. For fc = the initial conditions in the following sense 



t'o, 



A' 



A0,T, 



W^ 



WQ^- 



on fl. 



ipdx. (4.1c) 



(4.1d) 



In (|4.1dp . we denoted by Aq.t G i^'?(r2; R'^+^) and wo.t G £^(17) respectively suitable 
approximation of the original initial conditions Aq G L'?(17;M'^+^) and wq G L^(0) such that 

Ao,r^Ao strongly in L« (17; M^+1), and \\Xo,r\\L^-Hn-M^+^) < Cr-^^^'^'^'K (4.2a) 

wq^t -^ Wo strongly in L^ (SI), and wo.r <S ^^(f^)- (4-2b) 
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Moreover O^^f..^ G ^^(r) and 5^ G L°°(r) are defined in sucli a way that their piecewise constant 
interpolants 

[^cxt,r,&r]W:=(Ct,r,&r,) for (fc-l)T<<<fcr, k^l,...,Kr. 

satisfy 

^oxt,r ^ 6*0x1 strongfy in L^(S) and 6^ ^ 5 weakly* in L~(E). (4.3) 

To see why (|4.ip indeed forms a correct time-discretization of the weak formulation of 



(|3.1[) . note that (|4.1ap already contains the flow rule since, if we can find a minimizer, the first 
order optimality condition in A evaluated at Aj: yields the discrete version of (|3.3p . 

Realize also, that we have added the regularization term rjAp'^ to the time-incremental 

problem, cf. (|4.1ap . This shall assure that even | "^ ~^ j'? G L^(f2) (although it does not hold 

uniformly for r — >■ 0) and in turn also the existence of solutions of the enthalpy equation in the 
classical weak sense. Of course, this regularization term will be shown to converge to as we 
pass to the limit r — > 0. 

Proposition 4.1 (Existence of discrete solutions). Let p.Sp hold and let also w° > 0. Then there 
exists a discrete weak solution according to ()4.ip such that w^ > for all k = 1 . . .T/t. 

Proof. First note that (j4.1cp is decoupled from (j4.1ap and (j4.1bp . It is easy to see that, for 
any ly G ^^(51; M''), there exists a unique Um G H^{R'^) that solves (j4.1bp since to = id^i^ G 

L^(r2;IR''). Therefore we may proceed by the direct method to prove existence of (j4.1a|) . (j4.1b|) . 
i.e. take at time-step k a minimization sequence {qkj}j^o — {{>^k,j,^k,j)}fLo of (|4.1ap . Due to 
the coercivity of the cost functional (thanks to assumption p.5ap ) this minimizing sequence 
converges weakly* (at least in terms of a subsequence) in the space L^{n;M{R'^)) x L''{^; R'^"'"^) 
to some Qk- Moreover, note that, again due to assumption (|3.5ap . Vk is a Young measure. Then, 
by the convexity of the functional in A, v, by the fact that mkj -^ ruk in L^(i7;IR'') (because 
Vkj as well as Vk are in '3^P{Vt\ W^) and Vk.j -^ Vk in L'^{Vl\ M.iW'-))) and the boundedness from 
below of (j), Qk is the sought minimizer of ()4.1ap at time-step k. 

The existence of solutions to (|4.1cp for fc = 1 (and subsequently also for all other k) can 
be proved by standard methods exploiting theory of pseudomonotone operators (here rather 
trivially since the problem is semi- linear and thus the underlying operator H^{n) -^ iJ^(J7)* is 
weakly continuous). Note that the right-hand side can be represented as an element of H^{n)* 
due to the integrability of the initial data, the suitable choice of the time-discretization of the 
external heat flux ^^^t t G ^^(r), and thanks to the regularization term t|A|^'' in (|4.1ap which 
causes |^— ^^ |'J g L^(J7), as already mentioned above. 

Let us test (|4.1c|) by [w'^]~ = min(0,w^), which is an element of H^{il) as w'^ G H^{n) 
and hence a legal test function. We get 



\k \k—l . \k \k- 



'U] 



wlr]ds. 



(u;^[«;^r +r/C(AJ:,^«^)V«;^ V[^^]")dx- < ^^ (^5*(^^A^) 

Now Q(uj^)\w\\' = 0; here wc realize that we have defined Q(uj) = for w < 0. Further we 
realize that Jg C^""^''' ) > and g| ^""^"'' |g > q which implies that (5^( ^"~^"'' )[w^]~ < and 
e\—^—^ — [w^]^ < 0. Using these facts and by exploiting that b^Ocxt.Aw^]^ < (a consequence 
of ((33f| ). we get 

/ {[w'^n^ + TKo\^[w';]-\\\x < / u;^-i[w^-]-dx. 
Jn Jo, 

Using this equation recursively and when also taking into account that w° > gives w^ >0. D 
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Let us introduce the notion of piecewise ajfine interpolants A,- and w-r defined by 

t-{k-l)T ,,^ i,, kr-t 



[x.,w.]{t) := ' ^^ ''' {x'.M) + ^^—{x';-\w';-') 



for te [(fc-l)T,fcT] 



with k = 1,...,T/t. In addition define the backward piecewise constant interpolants Vr, Ar, 
and Wr by 



[Dr,~Xr,Wr]{t) := (l^^,A^,W^) 



for (fc-l)T<t<fcr, k^l,...,T/T. 



(4.4) 



Eventually, we will also need the "retarded" enthalpy and magnetization piecewise constant 
interpolant w^, m^ defined by 



[2£,(t),m,(<)]:=K"Md.z.^^] 



for {k-l)T <t<kT, k^l,...,T/T. 



(4.5) 



Note that from now on we use C as a generic constant, which ma change from expression 
to expression, and do not specify its dependence on the problem parameters such as e,q,p, |il|. 

Proposition 4.2 (A-priori estimates). Let the assumptions p.Sp hold and let t < tq for some 
To > fixed. Then the interpolants of discrete weak solutions satisfy 



sup J^\ ■ \^ tiyrdx < C, 
te[o,T] 



'^'^llL9(Q;R<i+i) - ^' 

X II < r'T-~i/29 

^llL°=([0,T];L29(O;Kd+i)) - ^^ ' 

'^^llL°=([0,T];Li(a)) - ^' 






with any I < r < 



d+2 
d+1' 



(4.6) 

(4.7) 
(4.8) 
(4.9) 

(4.10) 
(4.11) 



Let us emphasize that our strategy of proving ()4.7p and ()4.9p will slightly deviate from 
the standard approach based on testing the flow rule (j2.10bp by A together with the entropy 
equation (|2.6p by 1, which, when adding these two, would lead to canceling of the dissipative 

heat rate as well as the calorimetric term Q{w)a-X, see e.g. [34]. In contrast to p^ however, 
here we use the retarded enthalpy in the time-incremental minimization problem and hence 
the strategy of [34] would not work in our case, as the calorimetric term would not cancel out. 
Therefore we exploit the rate-dependent dissipation term that yields more regularity than the 
rate-independent contribution. 

In what follows, we will use the abbreviation ((•, •)) for the scalar product in H^^{il; M''+^); 
in view of the specific choice of the norm on iJ~^(ri;R'^+^) in Sect. [21 this means 



(X,v} := f VA-^X-\/A-\dx. 



(4.12) 



Proof of Proposition \4-S\ Before giving a rigorous proof we give a formal heuristic sketch, using 
the system (|2.13[) . on how estimates (|4.6p - (|4.1ip can be established; we shall always point to the 
adequate step in this proof where the formal procedure is performed rigorously on the discrete 
level. For this heuristics only, we shall assume that all functions are as smooth as needed. 
Furthermore, let us denote by (i^. A, 9) the solutions of (|2.13p . 

First of all we exploit (|2.13b ) , which can be reformulated as minimizing part with respect 
to i^ the magnetic of the Gibbs free energy © defined as 



e5(<,i/,A) 



(pmv — h{t)-mdx + 



-\Vurn\^ dx + -\\X-L .iy\\^_,^^^,^^^,y (4.13) 
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Then just using the chain rule and reahzing that the partial derivative of with respect to v 
evaluated at the minimizer has to be zero leads to 

©(to, i^(to), A(to)) - ©(to, i^(O), A(0))+ y" "^©(t, v(t), A(t))dt 

= <3{tQ,u{0),X{0))+ f {h{t)m+K{\-L.v,\}\dd, (4.14) 

for somec arbitrary tg e [0,7"]. Using the Young inequality and assumption p. 5b ) on the 
coercivity of the isothermal part of the anisotropy energy and also (|3.5b ) yields (cf. Step 1 
below) the following estimate 

C I l-r.i'(to) < r I (^|A|« + C|-|P.i^)dxdi + C. (4.15) 



Jn 



Inequality ()4.15p needs to be combined with an estimate on A. In order to get it, we 
multiply the flow- rule p.l3ab ) by -!-A and integrate over f2 and [0,to], which yields (note that 

r I -{5l{\) + -\\\'')dxAt^ /" (-((A-L.j/,A))-- f {Q{w)^e^)a-Xdx)dt. 
Jo Jn Q^ 1 ' Jo ^ Q Q Jn ' 

Using the chain-rule, we can rewrite this as 

/" / i(<S^(A) + -|Ar)dxdi+-||A(to)|lff-i(0:E-+i) 

< -\\m\\H-Hn:R^+^) - r' {-{L,u,\} + - f {e{w)~e,)d-\dx)dt, (4.16) 

which, by usage of the Young inequality (cf. Step 2 below) , yields the following estimate 

r / -(<5^(A)-f-|Ar)da;dt+-||A(to)||H-MaM-^^) < C \ (f |Ar+C|-r.z.-HC|«;|) d.Tdi+C. 

(4.17) 
Finally, multiplying the enthalpy equation ()2.13b ) by ^, integrating over D, and [0,to], a-nd 
using the Young inequality as already above gives 

-^ r I wdxdt< f" I (^|A|« + C|u>|)d.Tdi + C (4.18) 



Adding (|4TT5l) . (|4Tf| and (|4TT8l) gives, on this heuristic level, (j^G]) . (|4J| and (jlJ]) . 

Estimate (|4.10p is got by testing the enthalpy equation (|2.13b ) by 1 — ,^2^\a (cf. Step 4, 
below) while (j4.1ip is got from the enthalpy equation itself. 

For clarity, let us divide the formal part of the proof into five steps. 
Step 1: Using the time-incremental m,inimization problem,: In this step we perform the procedure 
that on the heuristic level led to (|4.15p : together with Step 2 and Step 3 it will give (|4.6p - ()4.9p . 

As we know that (i/^, Aj.) solves the time-incremental problem (j4.1ap . we may write 






\[-x[-' 



dx 



, , A^-A:, -\ re 



T ^T 



dx, 
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which can be rewritten using the magnetic part of the Gibbs free energy as 



<0(a-i)r,z.r\Ar') + 



(l-l)T 



e',{lT,,.i) + &'M,K))dt 



where the inequahty on the second hne is got by the discrete chain rule (relying on convexity). 
Realizing that ©((/r, i/l) = J^ hr{lT) ■ rnJ'^Ax [hr denotes the piece- wise linear approximation of 
/i)and(5'^(i/i,Ai) = >.K-L.z.[, 



/ A'. -A- 



and summing from to fc gives 



e5(ifc,z^,(ifc),A.(ife))<0(O,i>r(O),A,(O))+ / n hr-fhrdx + }<{K-L,Dr,\r})dt (4.19) 

with tfc = fcr, a discrete analogy of (|4.14p . Exploiting ones again the discrete chain rule as 

1 - 1 - /■*'=._ 

2l|Ar(0)||^-l(f2;Rd+i) - 2l|Ar(ife)|lH-i(0;R=i+i)da; > - / ((A^,A^))di, 



we can rewrite (|4.19p as 

<&{tk.Vr{tk)X{tk)) + '-^\\K{tk)fH-^(n-M^+l) <^l|Ar(0)||^-i(n;R<i+i) 

+ ©(0,Pr(0),A^(0))+ / ( / K ■ rfirdx - >f{L .Vr,\ri)(it 



Estimate on the right hand side (using Young inequality and (J33J:)) 



/i^-?TiT-da;dt 



Jn 

tk 



.p^)d.Tdi < C 1 + 



^L ti'r, AT-))di 



< C f"! {\K\^' + I 
Jo Jn 

Jo 



' • VrAxdt J 



Jn 



< 



< 



J V4gCem ^^ '^^'^''^'^ + ^ll^"^'r||^-l(f2.Rd+l)jdi 



\K\\l 



L<!(0;R''+1) 



c 



' • i^T-dx ]dt + C 



with p' = ^ and £„„ the specific constant for which WKW^-^^.j^d+i^ < Ccm\\K\\1^Q.^d+ 
Note that we also estimated, thanks to g > 2, ||i»i^r H^-if^.^d+iN < C(l + \\L»Dr\\1^2(Q.Ts^d+i) 
Cil + J^\L\^*Prdx)<C{l + J^\-\P.Drdx). 

On the other hand, we may estimate (3{tk, Vr{tk), Xritk)) from above by using p.Sb .'l as 



< 



&itk,Vr{tk),K{tk))> / (t)'iyr{tk)dx 



h{tk)-mr{tk)dx 



> 



>C 



ct{l + \-\'''Mtk))dx- I (^C\h\P' + ^\-\P.Dr)dx 

cU \-\P,Vr{tk)dx-l) (4.20) 



Combining (|4.20p with the estimates on the right-hand side of (|4.19p one gets 



C \-Y''Vr{tk) + - / |A,(ifc)p9dx< 

Jo ^ Jo 

a discrete analogy of (|4.15p . 



Jo 



\^\i + C\-\P»Vr)dxdt + C, (4.21) 



Step 2: Testing the flow-rule by X^^: In this step we perform the procedure that on the heuristic 
level lead to (|IT7|) : together with Step 1 and Step 3 it wih give (|TB))-(| 
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First note that, as (z^^, A^) is a minimizer of (|4.1ap . the partial sub-differential of the cost 
functional with respect to A has to be zero at A^. Realizing that this condition holds at each 
time level and summing these conditions up to some k leads to 

(d*s{Xr) + -\K\Adxdt < [[K{\r~L,Vr,Vr-\r} 

+ j ((eK)-0c)a-(wr-Ar) + 2<zT|A,p9-2A,(z;^-A,) + <5JK) + -|«.|9)da;')dt, (4.22) 

where Vt is an arbitrary test function such that Vri-^x) is piecewise constant on the intervals 
{tj-i,tj] and Vr{tj,-) G L^gj^^.^d+i) foj. g^gj-y j. As 

1 - 1 - /■*'=._ 

2ll^^(0)llH-i(n;R<'+i) ~ 2\\^r{tk)fH-i(n-s.'^+^)^^ >- (Xr,K}dt and 

/ |Ar(0)|^'' - / \~Xritk)\^''dx >-2qf'[ A^A^|A^p9-2dxdi, 
Jn Jn Jo Jn 

hold by the discrete chain rule (thanks to the convexity of the involved functions on the left- 
hand-side), we may test ()4.22p by Vr ~ 0, which effectively executes the test of the discrete 
version of the inclusion (|2.10bp by -X^-. This yields 

' I {5*s{K) + -|A.r)dx-dt + |||A,(ffe)||?,_,(f,^R.+,) + / r|A,(tfe)pMx 
Jo 9 ^ Jn 

< / f >^l|-^»i>r||ff-i(0;R<i+i)||Ar||ff-i(n;R''+i) + / (|6(w^) -6'cl |a| |A^|)da;jdi 

+ fllAr(0)||l,-,(O;K^+,)+^r|A,(0)p«dx 

+ f (cie(«;j-0er' + ^iA.r)dx')dt+c 

-/ y (^' ■ I" ""^ + ^'"'^1 + ^\Xr\'^)dxdt + TWXoAll.^n-M-^^i^ + ^' (4-23) 



by applying Young's inequality to the terms ><^||i»i'r ||_ff-i(0;R''+i)ll'^r||H-i(n;R''+i) as well as 



0(^t) ~ ^c||a||Ar|- Subsequently we estimated \\L*h'r\\''^-i,^.gd+i) the same way ay in Step 1; 



Com was also chosen like in Step 1. Eventually, \Q{w^)\'^ < C(l + \w^\'^ /") < C(l + \w^\) due 
to dHU) and q'/LO < 1 (cf. IKEH)). 

Step 3: Testing the heat equation by 1: In this step we perform the procedure that on the 
heuristic level lead to ()4.17p : together with Step 1 and Step 3 it will give (|4.6p - (|4.9p . Summing 
the discrete version of the enthalpy equation (|4.1cp from to tk leads to 

/ I / I Wr<y3 + /C(Ar, U'r)VWr 'Vl/Jldx + / 5r 6(Wr )</'d<S' J di 

= [ ( [ [si^iXr) + elKl"^ + e{wr)d-Xr^ipdx + f bJ^^t^^ifdSjdt, (4.24) 
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where ^ is an arbitrary test function, such that Lp{-^x) is piecewise constant on the intervals 
[tj-i^tj] and tf{tj, •) G H^{Vl) for every j. Then we test equation (|4.24p by 1 to get 



( I WrAx+ I hrQ{Wr)ds\dt < / ( / UHK) + e\\r\'^ + \Q{Wr)\\a\\\r\]dx 

+ j |6r^ext,r|d5'jdt. 

Estimate the third term on the right-hand side similarly as in Step 2 to arrive at the expression 

/ I Wr{t,x)dx+ I bre{wr)dSdt< j 2\Xr\'^ + C\wr\dxdt. 
Jo Jn Jr JQ 

Multiplying this by e/{8q) and adding to (|4?23| and (J4?2T|) already yields (i4J|) - (j4J| by the 
usage of the discrete Gronwall inequality. 

Step 4- Estimation of\Iwr'- In this step we prove (|4.10p . 

Let us test (|4.24p by ri{wr) where r]{w) = 1 — ^_/^^„ with a > 0, which, due to the non- 
negativity of the enthalpy is a legal test. Notice that due to the discrete chain rule relying on 
the convexity of fj 

/ fi{T) — 77(0)da; = / —fi{wr)dxdt < / Wr'ri{wr)dxdt, 
Jn Jq di Jq 

where fj denotes the primitive function of r] such that fj{Q) — 0. Realize also that due to the 
fact that ■q{'w) > also f^T) > and hence we may write 

hqg / — — -\i, dxdt =^ kq / \'ViVr\'^'r]'{wr)dxdt < / /C(Ar , WT)V'iDr-Vwi-??'(wr)da;di 
Jq (l+lf) +" Jq Jq 

/C(Ar, WT)'VwT-'Vr/{wT-)dxdt 



IQ 

< / fl{Wr{T))dx + j IC(Xr,U)r)'VWr-^ri{Wr)dxdt + j b.rQ{Wr)ri{Wr)dSdt 



&T^oxt,r??(wi-)d5dt + / fi{wo)dx + / {6g{Xr) + Q{wr)a ■ Xt + e\XT\'^)'r]{wr)dxdt 
E ' Jn Jq 

<C+\\f,\\mQ) (4.25) 

where we used the obvious bound \r](wT)\ < 1 and abbreviated 

fr := S*g{Xr) +e\Xr\'' +Q{'Wr)a- Xr. 

As to the L^ -bound of fr, realize that 

f S*siXr) + e\Xr\'^dxdt < 2e||A^||^,(Qj + C and 

Q{Wr)a- Xr dxdt < C\l + ||Wt(0IIl1(Q) +eII'^r|li,(Q) 



similarly as in Step 1 or Step 2 of this proof, and use the already shown estimates ()4.7p . (|4.9p . 
Thus, (|T^ yields /^ J^fyJ^dxdi bounded. Combining it with gH) like in [31|1], cf. also pi 
Formulae (4.29)-(4.33)], we obtain (|CTI)) . 
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Step 5: "Dual" estimate for the time derivative: Notice that 



N I. i._i N 



k i- — 1 



= > sup / vdx 

N 

= V sup / (-/C(A^,«;^■)Vw^V^; + (5J(A,(tfc))l. + e|A,(ife)|'^« 

^.^^«ewi.°°(n), ||u||<i Jn ^ 

+ Q{w'')a-\ritk)l^dx+ f {br9<,^t,r-br&iw''))vdS 
< sup / ( — IC{XT,WT-)VWr -Vv + 5g{Xr)v 

■flec([o,T],iyi-°°(o)),||t)||<i JQ ^ 

+ e\Xr\''vQ{iUT)a ■ XrVjdx + / (brOcxt.r — br'd{wT))vdSdt. 

Now because of all the preceding steps we may use the Holder inequality for all terms on the 
right-hand side to get estimate (|4.1ip . D 

Proposition 4.3 (Convergence and Existence). Provided ()3.5p holds, there exist {ly, X, w) e 
(^P(n;R''))[o^'^l X M^i'«([0,T];i«(fi,R'*+i)) x L''{[0,T],W^^''{n)) and a sequence t -^ such 
that 

Xr-^X in i°°([0,T];i«(fi,R'*+i)) and A^(i) ^ X{t) Vi G [0,T] m L'(fi;M'*+^), (4.26) 

A^^A mi«(Q;M''+i), (4.27) 

Wr-^w LmO,T];W^^''{n)), r < !^ and w^ -^ w in L^{Q), (4.28) 

fhr^m in L'^{Q;K.''), (4.29) 

L.Vr^L.v inL^{[Q,T]\H-^{VL-M.'^+'^)). (4.30) 
Moreover, for each t G [0,T] t/iere exists a subsequence Tf^m such that 

^r,„,(t)^KO *nL-(17;A^(M^)). (4.31) 

Every (i/, AjW) obtained in this way is then is a weak solution of (|2.10p . 

Remark 4.4. Note that the Young measure i/ obtained by (|4.3ip surely does not need to be 
measurable as a function of time. However exploiting the convexity of the magnetic part of the 
Gibbs free energy and, in particular, its strict convexity in the moments m = idtiy and Ltiy 
we shall prove that, in contrast to v, the moments m and Lmiy are measurable in time. 

Proof of Proposition \4.3\ Again for lucidity, we divide the proof into six steps. 

Step 1: Selection of subsequences: By the a-priori estimates proved in Proposition 14.21 we may 
find a sequence of r's and (A,^) G (Vi'''([0,T]; L'?(f7,M'*+i)) x L''{[0,T],W^^''{n)) such that 
(|4:26l) - (i08l) hold. Moreover, there exists S G L'^{Q;W^+^) and h G L(2[0,r]xK'';K'') such 
that 

L*Vr ^'B. in L^{Q- M^+i), (4.32) 

Vu^^^h in L2([0,T]xR'';M'^). (4.33) 

Indeed, by (j4.7p we know that A,- is (considering also the integrability of the initial condition, 
cf. p. 5^ )) bounded in VF^''([0,T]; L''(rj;M'*+-'^)) and hence converges weakly to some A in this 
space (for the subsequence selected). As W^^'i{P,Li{Sl,W^+^)) C C([0,T]; L9(17,M'^+i)) we also 
get that A^(t) ^ A(t) for all t G [0,T] in L9(0,]R''+i). 

Moreover, At converges weakly in L''{Q; R''+^) to A as well, because jjA^ — Ar ||Lg(Q;Md+i) < 

T||A||L,(Q;Rd+i) -^ 0. Further, again by dH?]), A^ is bounded in L°°{[0,T];Li{n;R'^+'^)) and 
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also BV([0, T]; L''{n; R'^+^)). Hence, by making use of a slight modification of Helly's theorem 
[231125] Xr{t) -^ X{t) for all t e [0,T] in L«(0;K'*+i), too. Note also, that due to the fact that 
q > 2 and the compact embedding L^^f^.^rf+i) ^ i/-i(f7;M'*+i), Kit) -^ X{t) strongly in 
i?-i(0;R'^+i) for aht G [0,T]. 

For this sequence of r's, L»Dt- is bounded in i^(Q;R'*+^) due to ()4.6p : therefore (|4.32p 

follows just by Banach's selection principle. Similarly, also rhr is bounded in L^(Q;M'^) due 
to ()4.6p and hence, by standard theory for the elliptic equation ()4.1bp . Vu™^ is bounded in 
i2([0,r] xR'^;M'^) ; this readily gives (|T33ll . Due to gS]), Wr is bounded in L'-([0,r]; W^'^'in)), 
^ < ^+1 3'iid therefore converges weakly to some w in this space. Having the dual esti- 
mate on the time derivative of Wr (|4.1ip . we exploit the Aubin-Lions- lemma generalized 
for measure- valued derivatives (see (22) to get that Wr converges even strongly to w in 
^(d+2)/(d+i)-/J([o,T]; yyi-/3^('i+2)/(rf+i)-/3(n)) and after interpolation L^~P{Q) for any ^ > 
small, so that the traces converge strongly in L(''+2)/(^+i)-/3([0,r]; L('*'+'*-2)/{d=^-2)-/3(p)). f^j. 
any small /? > cf. [Ml Formulae (4.42) and (4.55)]. Moreover the estimate (|4.1ip assures that 
Wt, cf. (|4.5p . converges strongly in L~d [Q) to the same limit as Wt- 

Thanks to the growth condition |6(iu)| < \w\^/'^ and assumption p.5ep we have that 
8(w)| < \w\'^/'^ < jwl^/* . Hence exploiting the continuity of Nemytskii mappings in Lebesgue 
spaces and using continuity of 6, we have that 6(u't-) -^ Q{w) in L* [Q). Similarly also 
<d{w^)^Q{w) inL'j'iQ). 

Now, take any i e [0,T] arbitrary but fixed. Then due to the bound (|4.6p select a subse- 
quence rfc(4) such that Vr^^^At) -^ v{t) in L^(r2; A^(IR'')). As indicated by the index k = fc(t), 
this selection may depend on time t. Estimate (|4.6p then also assures that z/ is a Young mea- 
sure. Since the growth of L is strictly smaller than p it holds also that L»Dt^ ^^{t) -^ L*v{t) 

in L'^{Vt;W^+^) and therefore also L»Vr^^^^{t) -^ L»v{t) strongly in H-^(yt;W^+^). Similarly 

Vitid < u^ (t) —^ Vwid . u{t) in i^(f^; K'^)- Note that this does not necessarily imply that L • i/ = S 



^,(i) -- V-«id.i.(t) 



or that Vuid mu ~ h (where S and h were defined above) . 

Step 2: Minimization principle: Let t be still fixed. A direct consequence of (j4.1ap is the discrete 
minimization principle that reads as 

for any z) G ^p(J7;R'^) with © defined in (|4.13p . Here we denoted tr^^,^s = I ■ Ti^{t), where 
I = minsgN{i < sTk{t)}- Applying liminf,-^ _^o on both sides and using the definition of the 
magnetic part of the Gibbs free energy we get that 

liminf®(^rfc„,i^r,(,)(i),Arfc„(t)) < 1™ ©(tr,(,),I>,Ar,„(t)) 

p.i^- ti[tr,„J-max+\\Ar,,,^(t)- L,iy\' 

where m = id • i> because of the continuity of h and the strong convergence of Xt (t) in 
i7^^(J7;R'^+^). As to the left-hand side, because of the boundedness from below of (j), we may 
estimate 



.^"^io / '^•'^"'*(^^Mt))'"^^-^+II^^M*)W"^"'L-l(0;E<*+l) = ^(^''^'^W)' 



{(jyv{t) - h{t)-m{t)) dx + K\\X{t) - L»v{t)\\\_^ 
n 



- - 11- _ ||2 

< liminf / 0.z^^,(^)(i) -/i(i^^(^j)-m^^(,)(i)da;4-ii:||Ar,„(t) -i.i/r,„(i)||ff-i(n-Rd+i), 

which already gives the sought minimization principle. 

Step 3: Measurahility of Lmu and Vitm(')j strong convergence of Ltv: Let t remain fixed as in 
Step 2. Then, the convexity of 0, the convexity of the space of Young measures and the strict 
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convexity of © in the terms A— L • v and Vm„i, used for m = \A»v with v corresponding to some 
minimizcr of p.ip and A = A(<) already selected in Step 1, ensures that both \—L*v and Vu„i 
are determined uniquely. Then, since A(i) is fixed, also L*u is determined uniquely, although 
the minimizer v does not need to be. 

In turn it means that L»Vr{t) -^ L*v{t) in L^(il;M^+^) and VUf^^^t) ^ V"m) in 
L'^{yi;W'-) for the sequence of r's already selected in Step 1. Hence, by usage of the Lebesgue 
dominated convergence theorem, L»i'{t) ~ S(t) and h — Vu„i(t) for a. a. t G [0,T]. This in 
particular shows that both L • i/ and Vu^u\ are measurable. 

Note that the above implies also \\L»Dr{t) — L»j/(t)||^-i(Q.]jd+i) -^ for all t G [0,T]. 
Realizing moreover that sup^grg j^i \\L*VT-{t) — L*iy(t)\\'^_i,^.g^a^is^dt is bounded by a constant 
due to (|4.6p we may use Lebesgue dominated convergence theorem to get (|4.30p . 

Step 4- Convergence of the flow rule: We arc now in the position to pass to the limit in the 
discrete flow rule, i.e. (|4.22p but integrated to T instead of t^. First choose a test function v e 
L''{Q;M.'^~^^) and consider its piecewise constant approximations Vr such that zv -^ v strongly 
on L«(Q;K''+i) and moreover ||wr||L2,(Q.Rd+i) < Cr-^/^^q+i)^ gy convexity of S*s{-) + || • |«, we 
get that 

'' ((5S(A) + -|A|«)da;di<liminf / {S*siXr) + -\Xr\'')dxdt + f T\\r{T)\'^''dx. 
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As to the convergence right-hand-side of (|4.22p . we use that Q{w^) -^ Q{w) in L"^ (Q) 
to pass to the limit in J„ 9(w^)— 6'c)a-(fr— Ar)dxdi and (|4.26p as well as (|4.30p to establish 

the convergence of /„ x(l\r—LtDr,Vr~Xr))dt. Since ||Ao,r||L29(f2;Rti+i) < Cr~-'^/('^'+-'^' the term 
T j^ |AoP''d.T converges to zero. Similarly 2(7r|Ar P'^^^Arfr can be pushed to zero thanks to (|4.8p 
and the blow-up for Vt specified above that allow us to estimate | /„ 2(jr|Ar P'^^^Arfrlclxdi < 

2gT||Ar||/2,(Q.Rd+i)|lwr||L29(Q;Kd+i) < Ct^i'^+^i. Altogether, applying liminf^^o to both sides of 
the discrete flow rule, we get 

f 5*s{\) + -\\\'^dxdt< f (>c{{X-L.iy,v-X} 
Jq 1 Jo \ 

+ f (^{Q{w)-e,)a-{v-X) + S*s{v) + -|w|9)dx')dt, (4.35) 

for any w £ L'i{Q;R'^+^). 

Step 5: Strong convergence of Xr'- First test the discrete flow rule (cf. ()4.22p . reformulated using 
the convexity of \ ■ |') by As^i- being a piecewise constant approximation of the function A such 
that As,r -> A strongly in L«(Q;]R^+i) and moreover |lAs,^||i2,(Q.RA/+i) < Cr-i/t^-J+i). We get 

[ 5*s{Xr)dxdt+ j T\Xr{T)\'^'^dx< f T\Xo\^'^dx+ f (x{{Xr-L.Dr,Xs,T-Xr} 

Jq Jn Jn Jo \ 

+ J {^Ukr) + e\Xrr^XriXs,r-Xr) + (e(w,) -^e) A' (Ag.r- A^ + 2qT\Xr\^'>-^Xr-Xs,r)dx\dt. 

(4.36) 
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Symmetrically we test the continuous flow rule reformulated as above by A,- to get 



5*s{\)dxdt< I (x{\~L*v,\r-\} 

+ / (e|A|'^"2A-(A^ - A) + {Q{w) - B^)d-{K - A) + (5J(f))dx jdi. (4.37) 
We add ()4.36p and (|4.37p . apply Holder inequality and limi-^o to estimate 

e lim(ll^r|li,(Q.ud+i) - II^IIL(Q;R'i+i)j(ll^^lli'(Q;R''+') ^ ll^llL'J(g;R''+i)j 

< lim e / (\\^\''-'^K - |A|«"^a) ■ (A^ - A)dxdi 

T-)-0 Jq\ J 

|A.r''A.(As,. - A) + 6*s{\s,r) - 5*s{\\ ) dxdi 




(I) 



((©(wZr)-^c)a-(Asr-Ar) + e(w;)a-(A^-A))da;dt+ / xllX-L.v^K-X} dt 
^- ^ '^ Ja- ^ ' 

(II) (III) 

+ / {x{\r-L,Vr,Xs.r-K}+ / 2gr | A^ p^-^^^ Ag.r dx) dt + /r|Aop« I <0. 



(IV) (V) (VI) 



When passing to the limit on the right-hand-side use that As,r — t- A in L^(Q;M''+i) to limit 
Term (I) to 0. For Term (II) we use the convergences established in Step 1, and the convergence 
of e(w^) -^ Q{w) in Li'{Q;W^+^) to see its limit being 0. When turning to Term (III), ^^HFJ) 
needs to be applied to find this term approaching again 0. Term (IV) converges to by applying 

(|426l) and (|4?30l) combined with As,r -> A in L'(Q; M''+^). Term (V) can be pushed to similarly 
as when converging the flow rule using the available blow-up conditions, Term (VI) converges 
also to by exploiting (|4.2ap . Passing then to the limit and using all above said, we arrive 

at ||Ai-||l<!(Q;e<*+i) ~^ II'^IIl<!(Q;R''+i)- Hence, by the local convexity of L'^{Q;W^^'^), the already 

proved weak convergence A,- —J' A in L''(Q;R''"*'i) turns to be strong. 

Step 6: Convergence of the enthalpy equation: Performing discrete by parts integration in (|4.Icp ) 
yields 



Wr"^ + /C(At, WT)Vu'i-V(y9)da; -f- / brQ{wr)'^dSAt — I l(Jo,T<y5(0)dx 

Q Jv Jo. 

+ / {5*s{\T)'f + e-\\r\'^^)AxAt + / Q{wr)a-\r'fAx + / brOcy^t.T^dSAt, (4.38) 

Jq Jq Jr 

for all (f piecewise constant on the intervals {tj-i,tj] such that ifi{tj, •) G H^{n), (p{T) = and 
If piecewise linear on the intervals (ij_i, i^], such that (p{tj, ■) = ip{tj, ■). Note that by such test 
functions we may approximate (strongly in the norm of U'(Q), p € [1, oo]) any ip G C^{Q). 

To make a limit passage for r ^- in this equation, we make use of (|4.28|) (and the 
approximation of ip mentioned above) to handle the term {UtP- Then use that At- — ;■ A 
strongly in L''{Q;R'^^^), which together with p.5ep gives IC{Xr,iVr) — > IC{X,w) strongly 
in any Lebesgue space, except for L°°{Q;M.'''^'^). Exploit also that Wr ^ w strongly in 
^(d+2)/(d-i-i)-/3QQ^y].p^i-/3,(d+2)/(d+i)-^(j^')) for any /3 > small (as shown in Step 1), so 
that the traces converge strongly in L(d-H2)/(d-Hi)-/3QQ^y]. 2^(d=+d-2)/(d=-2)-^(P)) Combining 
that with (|3.6p allows us to handle the left-hand-side boundary term. For the right-hand-side 
we exploit the strong convergence Ar — > A in L'^{Q]W^^^) to the limit in the terms expressing 
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dissipated heat. For the term Q{wT)a-Xrf> exploit Q{wr) -^ Q{w) in i^ (Q), for the right-hand- 
side boundary term we have weak convergence of 6r^cxt.T ^^ bOcxt in i^(S), which is enough to 
estabhsh the hmit of this term. D 



5. Less dissipative modification of the model 

In some cases it might be advantageous to generahze the presented model in the following 
way: assume that only one part of the internal parameter A, which we denote AX, to be sub- 
jected to rate-dependent dissipation as before while the other part A — AX evolves purely rate- 
independently; here A : R'^+^ — >• R'^+^ is a linear projection. Assume moreover that the coupling 
between the magnetic and thermic part is realized only through the rate- dependent part of the 
vector of volume fractions, i.e. Ker A C Ker(a-). We alter naturally the dissipation potential 

R(X) ~ [ 6*s{X) + -\AX\'^dx. (5.1) 

Jn 1 

We further suppose that the rate-independent part can be split as Sg{X) ~ Sg^ {AX) + St; (^^A) 

with ^-""A := A — ^A, and leave the Gibbs free energy unchanged as in (j2.4|) . The evolution of 
the system is then again governed by ()2.10p . Now, however, the system hosts a rate-independent 
process in evolution of the component A^A. In the spirit of j34j we can formulate this system 
weakly by exploiting the concept of semi-stability combined with energy balance: 

Definition 5.1. We call the triple iv,X,w) G (^P(f7;M^))[0''^l x BY{[0,T]; L^{n;R'^+^)) x 
L\[0,T];W^'^{fl)) such that m = id . j/ G L^{Q;R'^), L.v e L'^{Q;R'^+^) and AX e 

L«(r2;M''+i) a weak solution to (^1^ with e\X\'^~'^X in (|2.13bp and e|A|9 in (|2.13cp replaced 

respectively by e\AX\''~'^AX and e|AA|'' if it satisfies: 

1. The semistability; For all v in ^P(fJ;E'*), all X e H-'^{n;W^+'^) such that AX = Q and 
all t e [0, T] 

'^{t,u,A^X,Q{w)) <'^{t,D,X,Q{w)) + f S*s^{A^{X-X))dx. (5.2) 

Jn 



2. The rate dependent flow rule. (|3.3l) with AX instead of X and Av instead of v. 

3. The total energy balance 

©(T, u{T),X{T)) + f be{w) dSdt < (8(0, z^o, Aq) + / wodx 
Js Jn 

^cxt d5dt + / ©;(t, iy^,Xr) dt. (5.3) 

'E Jo 

4. The enthalpy equation.- p.4p with jQipT-L{dxdt) 4- J„ dg^{AX)(pdxdt instead of 

lo <5s(A)'/'da;di and J„ e|74A|^(y9dxdi instead of J„ e\X\'^ipdxdt; here we denoted % G M.{Q) 
the measure (—heat production rate by rate-independent dissipation) defined by prescribing 
its values on every closed set A = [^1,^2] x B, where B G fl is a Borel set as 

7{(A)=Var5j^(A|B;ti,t2). 

5. remaining initial conditions and the reduced Maxwell system p.2p . 

Theorem 5.2. Let p.5p hold and A be a projector R'*+^ — > R'*+^ . Then at least one weak solution 
to p.lOp in the enthalpy formulation in accord with Definition \5.1\ does exist and also p.7p holds 
but weakened as ^ G L^{I; W^'^{n)*) + M{Q). 
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Sketch of the proof. We proceed in a similar way as above by discretizing the problem in time 
and introducing a time- incremental problem (j4.1ap with the discrete version of |A|'' replaced by 

\AX\'- 

The existence of approximate solutions and also a-priori estimates are proved similarly as 
above with obvious modifications where there are needed; in particular Step 1 of Proposition 
121 gives the boundedness of A^ in L°°([0,T]; iJ-i(r2;R''+i)) n BV([0, T]; L1(0;M'^+1)). 

As to convergence, even in this case, we can establish obvious modifications of convergences 
(|4.26p - (|4.3H) (wc replace A by ^A when necessary). To obtain the semistability, we first note that 
a direct consequence of the time-incremental problem is the discrete semi-stability condition 



®(ir„„,i^..,„W,^'-A.,,,,(t)) < ©(i.,,„,i>,A)+ / <5^(A^(A-A.,<J)dx, (5.4) 



for any v € ^P(f7;M'') and any A € ii^-i(0;K'^+i) such that AX = 0. This semistability can 
be converged similarly as above, however as we do not know that Xrf.,^-. (t) -^ A(t) strongly 
in L^(S1;M'*+^) (would be necessary to establish the convergence of the right-hand-side term 
dg{A (X — At-^.j,)) we use that & is quadratic in A and employ the so-called binomic trick as 
in e.g. [Ml. This trick is based on choosing the test function as A = A — X{t) + Xr^^^^ (t) (where 
A is arbitrary but such that ^A = 0), subtracting the right from the left hand side of the 
semistability, use the binomic formula and then passing to the limit. 

The flow rule can be converged in the same manner as above, also we obtain the strong 
convergence of AXr — s- ^A in L'^{Q; W^^^). The only delicate part in this case is the convergence 
of the heat equation, in particular the convergence of the right-hand-side rate-independent 
dissipation terms. To be able to find the limit of this term, it is necessary to prove the discrete 
energy inequality 



(3{T,DriT),XriT))+ / T\XriT)\^''dx+ / {5*siXr) + e\AXr\'')dxdt 
Jn JQ 

<&{0,VriQ)Xm+ I T\Xrfl\"^Ax+ f <&[{t , V^ , Xr) ■ (5-5) 

Jn Jo 



To get this inequality, we exploit the convexity of the problem and realize that solutions of 
the original time-incremental problem are also solutions of the following auxiliary minimization 
problem 



Minimize ^(tfe, i/. A, e(u;^-i)) + / [t\X\^'^ + t5*s( ^—\ 



AK-AK 



9-2 



'A-A 



fe-i, \ 



AK-AXt 



AX-Xt"- 



subjectto {v,X) G ^P{n-W^) x L^'i{n-W^+^) 



dx 



(5.6) 



and vice versa. Inequality (|5.5p is then got by testing (|5.6p by {v^ ^,A^ ^). Using (|5.5p and 
also the "inverse limit energy inequality" 



e5(T,i/(r),A(T))+Var5j (^^A;[0,r])+ / (5*^(AA) + e|AA|«da;di 

> 6(0, 1^(0), A(0)) + / &\{t,u{t),X{t))dt (5.7) 
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that is a consequence of the semistabihty (see e.g. [34]) we can estimate 

/ n{dxdt) = Vai-s' {A^X, [0, T]) < Hminf / (St, {Xr)dxdt < hmsup / Jt, {Xr)dxdt 
Jq ^ ^^0 Jq r^o Jq 

<hmsup / (&'t{t,Dr{t))- f {S*s^{AXr) + e\AXr\'')dx]dt - &{T,Dr{T)X{T)) 
T^oo Jo \ Jn J 

< y ("©Jlt, z/(<)) - y {pi^ {AX) + e|AAr)dx') dt - ©(T, j.(r), A(r)) + 6(0, v^, Aq) 

<Var,j^(A^A,[0,T]), 

which yields Xvxir^Q \Qb*gJyXr)dxdt = Var^-. (A^A, [0,r]) and hence allows us to perform the 
limit passage in the heat equation similarly as above. Note that, due to the form of ©. &\ does 
not depend on A. D 
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